We estimate the age of a site by calculating the years since the last fire. We then fit a curve to model the recovery of vegetation (measured using NDVI) as a function of it’s age. An additional level models the parameters of the negative exponential curve as a function of environmental variables. This means that sites with similar environmental conditions should have similar recovery curves.
The model was last fit on 2022-03-17 13:19:55.
This version of the model was fit with 1076106 pixels including data from 2002-04-23 to 2019-12-19.
This repository was developed using the Targets framework.
── Attaching packages ─────────────────────────────────────── tidyverse 1.3.1 ── ✔ ggplot2 3.3.5 ✔ purrr 0.3.4 ✔ tibble 3.1.6 ✔ dplyr 1.0.8 ✔ tidyr 1.2.0 ✔ stringr 1.4.0 ✔ readr 2.1.2 ✔ forcats 0.5.1 ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ── ✖ dplyr::filter() masks stats::filter() ✖ dplyr::lag() masks stats::lag()
Attaching package: ‘arrow’
The following object is masked from ‘package:utils’:
timestamp
CmdStan path set to: /home/rstudio/.cmdstanr/cmdstan-2.29.1
These parameters represent the relationship of the following environmental variables to the recovery trajectory.
The plot below illustrates some example recovery trajectories. It currently just shows the top 20 cells with the most observations.
Compare estimated vs observed values for all pixels. This is not true validation - these pixels were included in the model fitting.
model_prediction %>%
# filter(cellID%in%cells_with_long_records$cellID) %>%
ggplot(aes(x=median,y=y_obs)) +
geom_hex(bins=50)+
geom_smooth(method = "lm",col="red")+
geom_abline(color="blue")+
# facet_wrap(~cellID) +
scale_fill_viridis_c()+
labs(x="Estimated NDVI",y="Observed NDVI",
caption = "Blue line is 1:1, Red line is least squares regression. Count is the number of pixels in that location",
title = "Estimated vs. Observed NDVI"
) +
theme_bw()
## Warning: Computation failed in `stat_binhex()`:
## `geom_smooth()` using formula 'y ~ x'
## Spatial Predictions
Maps of spatial parameters in the model.
```{=html}
<div id="htmlwidget-ba5659cc4d7ad7328234" style="width:672px;height:480px;" class="leaflet html-widget"></div>
<script type="application/json" data-for="htmlwidget-ba5659cc4d7ad7328234">{"x":{"options":{"crs":{"crsClass":"L.CRS.EPSG3857","code":null,"proj4def":null,"projectedBounds":null,"options":{}}},"setView":[[-33.998707,18.577],10,[]],"calls":[{"method":"addProviderTiles","args":["Esri.WorldImagery",null,null,{"errorTileUrl":"","noWrap":false,"detectRetina":false}]},{"method":"addRasterImage","args":["data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAC0AAABfCAYAAAB88E5QAAAIOklEQVR4nO2beZAU1RnAf+91z/TMzszeu+OyhGNBd7PgeiDBKyjxiIVGhVUMmFh4g5qKMaUmf5gyZQWrxKQwFY3RlGiquFRWxSMeiEYTDRgRCAFFsoC4suuy7MzOsTN9vJc/KAmJWhVnegIx8/tvprve96vvfdOv+3s9UKbM5yJKMahOLNXkM+iBXmT7bb7HkH4POGf2Un3w50l3vqw/79xCMf0c7OIrH9NkHS78SQVm3iJb2USTl/EzBOCjdOe1KzWuwokGUVKgpaCyP4sy/K9AX8pj9e5HNOz/gQilERo8UyK0BgFzZi/ztUSKltZ77jsgJJQmkHEAsMMmRs4jNJhn6bLZvqa7qPLILThfA7TW5HCtKrQUBBN5bMugtifFkhVzSnJ1KkramnsOL53+NCef4xCdNp/UqybZxgqev+u8ksh+QsHSaufC/VmeINn840vxbnVZ9UhnSWU/oeCalmNuFjo5wOiuZ8WdC2rQ4r/iW6ZMmTJlypQpU6ZMmTJfEPX+At+7R/8JRbcQ1AvXHRLxotCZJ7R66Yb/PfEvBfa9F2vvqStLln3fn/u9lZdrrzeNMbISEQmjU1nktDMQ1Zf6Fsv3/rTRuViYE0dC3oVw6MD3emiF7h9erLV+pegZ8F0agEgYOfUkcD1ERysIA4aTBw6rP/ygKPGSSMsTFgj1+puQySKaWrEjMXQuhWVEcFUeqqJFje/rTsDBiBENqO5eZPdGgi3Hot7eTEV6PeZlS4qu7ZJJA4hoELVxO3L3HsTx7ciqJmBJ8eMWr/b5qK13aOf5d5DVFjJmYVz0cLlLWaZMmS8Le7IPaPfhOSW/ty5qGXfVS1rrNVqt/p4evv3cA7IXLH6mpOIFSw/Zj35KrCE0BtnRwoUnJbj4qsdLJl6QdMreLzRk95FxExCvQ4RNxFPL4bhpvJsQfNwc89f0IApaVtWG23S2/VgUimjvbnZev4Yxt01k+LgTGXaT1Ml6hHVuyZbsL5xpteUODRBauhJDmAzFm6h6bBZi/AQqbMX6vYBp+S56MF9YWv91O8QiePuGyXlpLCMKSD4MZnl5MMVZI5pZ3t1bAtV/UtAUJuefqUMtUYwRUYxzp+NVNpB2BujJpGmPNmMbGsuYfviUB0DVr1eL/rds+lbtQ/9tPYaCkBGl1vJ4Y6CXdX0f++35LxT8EDDi++2I1gkMV9UCWQSSSKCW8VVJknZp15eiptBeNFNrpQmcczwiPgoCIfSmdchTf17Sm/2iVsTgjV0i0BZHNDSz5usrIJtAjD+ShRuWHJ4r4ieoZAa3rplJ71xG5mePMv2pUQC8O/jg4bUiHswzN/fxp457qNq7l/zCuQQtxfH1Ls2RGj/8PpOipadMNdi1LY9a+za7UjY3nb2XWssl4w764feZFC0dv3wc3+kah7rgEgJS80HaJGkbhIwYSq/RuQXna/XcvMPrvTx51q+EqIxg3/JzauZ30RJzOK26kvCiZbC668B5M65/wjdxXy5Ntve8DngalAtODu/Zp/cfcDxkTQS+2sKMRW1YGZtHH5pVdExfenmB3DBDpNDb/kI2ZJKZMR3n3QFy63oZXrMDxk6GSSZCwezvLi86475kWr3+Q01DDSJWDUe0ktc5zMcfo39VP/EbJiKmnE3KSxB9521Ey1GkqyqJBS8qOLYvmVbb96C3/h2n61X0y09iKQO3O8Hvl+Xxtu6GRA/RQB35SaeQqa4GIG2vLDjjvjQgRXUIb1cClXFwNvXApgcw6sJ0zssi6yqg8ghclcf2snjawdUOhjAKjueLtDHjIaHW3qrFUUfjLHkSlbI5OXk1zUaKxPowr0ztx6huIhKoIeXs5ZqbazDzXsHxfGuqi5ZWsGLkrplF+sZOavfl0VIwcusAZ/yiHikMtFZYMsIHzTGckMHM+V0FlYh//elQFN3/PrFoI7l7VnDdtbNpqnAYHdPc9EKWmdcP41omNd1J2gMpnGgQL1hYznzLtE724o2cSDoWJnTrVXwrtZ3J2qExNJb7prtMmDuME5Bk4xGcaBDpKqxkvqBYvt33Og9+W8t4DHFMK+z4EOfUs5DCwEztQ/ftQFTWgpC86IT45aJaDFex6p4LCorvX3mYAnvDHgLJLHLiWKQ0SeZ7qQvXod/bgRraijy5gzM3rqXt9imMjl17aK/TAEa8CmvKKETUIvPbP2NsfoNqqwlvxTK8niQYAnb1IMIWX0nk0LnCW2e+ZVpOv194T1+tdcbGdQSieRwCgZg2GalcRLwV8mn0hrWwaQu67qOCY/m6u+Vu7UcNO1iNJjrVh6yogsYj0cEw3vLFyNZRiBPPRL/3FrLjp4e+PACsW54URm0Ya/IIyNsgA7imgaNyyEvmIo49A4TBa3Vj6ZxX+DLu+46t2RYHT6G37UQ99zhmah9BgggnB33bsBf9htMsCKadgmOU5FFfvXmLJhREd/fgbOnHiEfAlOQ39vNxd4gxP2pDnnJ3wbFL1p9w1YtabluHu2Yj6fVJKsaFwFX7W2lfa4O8jTxhwSG+Tv8b8oUnoH08RjxC1XmViDHNiIZG9M6drAiNZ6iIwizd3viIBuzfvUKgrQ553kzu2pLgjwujeOZo0tEgr93+zYJnuXTSHacTOHoq927uY8ndFiG7niMSScycS6CxoqihS/OSCiDFN4QUBjfUDnD/NbswlCIfDTI4MkaiJlzc2D45fiZCTBMiFOWYzD7qP0oTzNrkLRNPFvf7L6k0gKi/Qoh4K+kpEazBPPUfDlE3kC1qzJJLA4jw+eLojkGcWJBA1iXa5///bktG57yVetYVn95/LFPm/51/AFKYGcGgIviuAAAAAElFTkSuQmCC",[[-33.941641403099,18.3003829808018],[-34.3371886093941,18.5267584324]],1,null,null,"alpha"]},{"method":"addRasterImage","args":["data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAC0AAABfCAYAAAB88E5QAAAHh0lEQVR4nO2bW2wcZxXHf983szve9XrXXntjO04anMglSoIIKSJtmiCKSAr0pgopAkGAoj6ABDwhhKoqV5IgghCIIpUHHqqAkKpWBJJyEaRVUhFBlNLSOJc2F9upL3G83l3be5vb9/EwcQhtitTd2SRF+3uxd2Sf8/eZMzPnMoYmTd4V0Qij6vj3NJ4HloW8a3foPmTYBrOVZ/T1n9XwPv1uP1srZpjGxktPawB7zTry9jiZln5KXj5MF0CIog+PPqOhgqts/jo6gxBx4sZl5tzQT2Y46aHe2K0B/nQpwaGROLOOIGEq8o7ByZxk3a6/hJoidUdandyuAWKmwtMGbRFY3m4zNGfxzJF22g/kOPa7jaFejHUZU0M/1BTLkM3zQv9d9LW6nM5bJKMay1Dcv/grDbk71RfpuRIi0wPdizh01sAutfD46hLrex9riNh5ajauBndolILOLv7gpjgyYbDvni82VOw8NV+IctU2gWkg+74lXhw3GEiFfjtu0qRJkyZNmjRp0qRJkxBQp3bcknal7rmHOrXz/ddnqVM7tDq39/0n/P8CNbxPq3N7Ghb90OcU6s09mmgEXBeUBiGgWEKuCW9O3Zih+uB2jWWB74OUUCwF4ienoacLuXpXXX7Dn8NCIDCbB6XA9YKIm/+ZwKlLP64rdRoiWq7cJljYDW1JMI1AeNUGKwqeD8qvy36om4D/wvdBGsH3tgM6Am2tyI/uqTslGycawLWDFIlG0BNTqHw1FLMNnXKq8z/QGBJmi1CuIO/Zd1Omqk2aNGlyE6juelDP7xUbSV0PF2//l7RIWog7eqE9ifOro8Hx32zR5hf2N+yeXHPtUd378DsiavS0QjaPbIvh/vLzDYt4TaLdpzdrAGdwCn+4ABUbhMS4bzXu8WFmDw4h4pFwlV5HTelh3LsCeXoIla+gPQWpBGRzMDmN0dNKcm0/8hM/uX3SQ72+TQOIDw8gkhb2mQKvfPUE5GfBijLy7Cy0xsJXeh3vPdJSBmVnxMQbmUFG4Py5Com9l7gy5bHhFys58MgrDZB6nYT3/AurtgmsKJQqRD+1EqOzhU2fkzg2uJ5GXZjgge8kGqH1GrUv9M/v1SgNI+OoywXyh6d49ZhEA+1Jg7UnDt0+OX2NqRxcfAs1kUfN2cQzmnhMMj3pkogbIUp8J/W9pHJ6p8bzIZtHF8vosoPsTCI3PtXQYr+uxlau2CoQAuIxxKoBZLoNopEbPnjCpP5uvFxBZ/M4zx6D9jbE2k8CcO7TDzVMeP2NbXcXoiuNOV3CPXqWSCbNzGmHv7/shiDvxtQfaduBSBSZacOfqjD0zaO89IKHW25chtQvumJDNguGxEi34JQlrq3pWGgy/bVNWr2+Teuxn4X6F4RylauT2zXJNvA9/MOvYixdAN3pYHaXSSM6u6AliegM55W3cMZicyUolaDqYHzsTsh0QCQCHSm4kkO/eR4iLejygVAiHoro6sGTMDIOaLicRQ+NBrneGoPOFIVfn4HKLBgRdOm3dQsP5XSNbv6sLs8YLLnfJLJ+OeQKkGqDdCo4C4VZdKGI6FsAjguOg9z485p9hxLpBRviZO6wmRssU95/HFrjQTUoRBDteAzRnoDCHJSDeZ568ds1RzyUAaR59wCplSWwXSb3DxNfsgjGJoKo+n4wVNeAKcH10IUKWLV3NqHVCOrs9zW+D5aFf+Q1zC1fBtOCyiz6ny9DqQpREwyJ3vAA0i4jYo/U5D+8oXrVBttFD55Dl1z0lfPos38Draje/XFYmAnSxfWQs1PgVNDF52pKkdDm03p4HNHfh0inMBf3gFb4/zjDqcdP8KEnehDr1sOdHeQokKuW+EBiMWZ5riZf4Q3VlcL547/AU0RWdCGyecZ/P0uxDM8/NsRnNl8k8dPvYulWLhUVL41nWZqsbY0RXk6f3K5x3GCAHo2gLkzgXiyAEMgWA+34RO/7IPRmyPf1cXQcHl26pSb/4UXaikLEDBZDuRlU0aE67pB8uB+R6UBfmsT+8xmsjZqOi6M8WkejEJ5oQSDYM2BBGtP3aUtEEcsWg5QIz8e4UsI7NYpzZprKzgd1bGttfWRodw858ITAV0G0LQs9U4RoEHXys/hvjDFysIo3Nkd0oAMha8/McLdben6tXA6+Ko2eygdl66olLOvO413IgtK0PHmwZtXh78bnHzIzRfSVaVSugmyNIDLtYBjgKxjox+5dRsx86BZfiPMIEayUM2lEPIZhTaF9/+pxA3q7oFqhpZCt3UWIcq8RRFsFIqUM0kapoA4RAip2Xf8915jd+PInhb741tWKToNx1Y0QwbGp3K2v8m6EWNwLpXIQbccNBANYUcTae8Eu1my7Ma9OAOq187CgM6jqxiYhV0B0d0MqBa6N6Pn6bfBweRuVE5OUD+cwTE16+ybQCmIpdHsfwnPqst2wSLc9dVgkB0zS3/gIjE4AAjoWIV2bIqW6bDdMNEBs6yFBNArJBKIlAXNT6IiFpr7etqGiAeTqnQLXBbMFfBfGB2mbK9dnMyRt/9vJmt0CKfFSGUSsHZ0buxluw0EXn9Nq+EfNtyWbNHkb/wbXBgWo9H39UAAAAABJRU5ErkJggg==",[[-33.941641403099,18.3003829808018],[-34.3371886093941,18.5267584324]],1,null,null,"gamma"]},{"method":"addRasterImage","args":["data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAC0AAABfCAYAAAB88E5QAAAIP0lEQVR4nO2ae4xU1R3HP+e+Zubu7HNm37vsAj5YhCCKkgJpSCuKWBuBEktQWqPBVzFWG9PGtNWmVpumsUkTglGx4m5rUVALViOlFotFNIJUcKUCurLue3d2dndm7sx9nP6xlZBYbZm5t0ozn/9m5s75ffO7v/s753zPhSJFPhURxKAPv9Muu1OC2ojkllnX+B5D8XvAu/d1yFM/X3Hndvlp1+aLr6LX75kU/FqPxpEPTDZtr+fclWN+hgB8LI+fH+iQ3SlIWyoLGx0W1WeoChlIJNWRb/taIr4MtvFwu0zmBFFdMidm01Kaoyel01hiU6pXUNr/EcqUO3wTXnB5eH+5XQK4EkIK9GV0njwa5axynUbLoPRPf/ZVMBQo2vvHzySAqcGX63PseqeEvf0CT8LOboGIXyeUKx7yvXtohfz5Wa2Z5bPqcao3sbBnFQ1LBjmSDHP5lLWBtNKPyXvw7R88IW0PlrfUw1g/KAqifHWgYj8m7/K4svVaEdEkQr1EIF3kWJ+fuooUKVKkSJEiRYoUKVLEJ+55o8N39+i/oWAL4advfj7CC+LHb3TIRzrbzzzh/xds6myXv3k3uOz77lM82tkuBaAIcCSoAvYPCX69aI1vsQIxVx7pbJeGAhkHDHVSdDIR4qOuKLouefHGywuK67upDpOZGLEEVWFIZgWaAlrYOfn76qe3FVQ6gYi+vu0a0RCV2B4IAbURuKDOoaIih5TguoWFLciA/CwmcoCAiCYZtwUzKl0aFoyxYlrhZzCBiYZJv7pUhzrTAyB2SokUQqAu58bD7TKqT3YQTYFV0/0/6SpSpEiRzwk5/Lj09t4V+Nq6oMnl6i3PyG/Om+CqKg052n3y+4yzXUa0KwPryXkvAlZf8+QnMyo9iJqECXFk9OHAMp6X6FU3PD0p6AWbHcd1vLJq6OpBnjhO/7R67j0wyO6eiK9CTyUv0eYKj0xFiFxEZ2qli/jwAGLGTEilqVPquPk8i3Uzg5uyT1v02h1PSYDqlRmSVWF2vRpnU6YeaadBCOTgMX5xMOy/0lM47QdxZkOOw90hIobE0RTuX9WDIuDvWpTWeRdS7pUwsH40CK0nyesW/rFrs+zPKORc2HvCYFadzYqpKQxVUm7ECatRdPWyL055ACxrWSvaKmxsD1KJENt3x9n6fgmGqhNWoygi0GV6/n36oYMmIz0lnOiPELZdtrxUy/vj/dx63gnOLm/0U+MnKOgW3r2vQ7ZVSi5ttrn/gMEHR8uY0zbGPRf5Zxf8OwraYd43f414fVBw244KLmtyCZfYHOwsI5nbEuhUXvBuvMmErnej/Oj3DVx9QYoHl/cA0JkIbkYs+Dam7GfluD3AvW+W0NVZycVzE9SbkpwnWD87mDIp+DE3HQVTa2bDojAT84exPcn+IZdDw8F1kIJHHiZBLBcFzyWasUDVMRSHAx/prLxpq1xwk0V3Gh5c4F/WC67peORbYvtgmnX7kiQigks2l9NammXxtBwzrrNOXvfCh5t9q3FfbLEHttVxcE8td+7WaDt/mGROpT8tGMqCrsAPzs8ybitsOOSP/euL6HPaEpRXODRVTjpII1mNeETSaIIiJPHeAeIhh9604IEDhR93+FJne/sekzu7DbpTYI0b1MRyLKzzSNmCV3oFw91RBgYi7LrN48hoN7dsnsJfv3tp3rF9ES1T22TKAF2EeK4rwW/3lDP17DF+eAEMWWnW/66J8bSG6wkMbdLXs13Ba99bkld8X/qSHDiK6bggBJc0NdG6dBBVQG/aoyaiY+VUlAmJp4ItFCKGS1VlLu94vtS0MvUugWEg6s4i66a4+bEWDicMTM3DUCJsfeuXCCkxohJz1GbJoiGuX5zIO55vvXMos1nGvCg7BlOcU2HRPWEwO+ZSaTSQckZ5e9jijvYmcOD527vw8OhPK8yJX3/aGnwTvfPE49L2BIdGVOJhSZkB7yUFlgtLm3M4EtoqFDzpMmZPdhlNSFrLbjxtDb4dX+zu1TA1jytb02RdWNzg8f05NcyslGRdQaNpM2TlEELQVFLPdFlCi6XmFcu3TH/t0edlNqviuIK6hhRjr2qUfsll9bwJxmzB5c0ulc+/jKiPMzi3jRq9EaEvzSu+b5lWFElZWQ4z7DDQa5JqDbF0dgpdkcyLZ6jq64cV19Iz5yyq97+Tt2BfRSMhndKI1adpmjpObW2GbX+r4L2kSlO0EhFrYXfvCY6PhRm+cDZyNP+Z0df17oon/iDjLePUhmHnSzXEEhaZVgNd85g1Z4SfXGTSkxrild4I9abHspb8XsP3ddFbXp3BzSkMCY/Wcyfo6zGRQ5LIhMPQFIP1L+e4uDFCLEzegiGA0621O56S2bRGLqsyPBwmk1YImx5h3cNxBVVVWZbMHWN+TZa51affoyGAc0TbUgFBKOxSU5MhndJpmJ7EyaqMD4dBSHa+VcbAueN5xwhkD7fmma1S0SSq7qLqHmbII6JBwhI4lkrfh6W8sG7ZF6B7nELH8pViYtRA+dfoliNwPBCKRCiSr16c4NY9+XePQEQDbFjRi+sIVN1DVSUhBcp0mB5zKTMkc2P57wUC2zLXpRQem9IJUtI3rY5fvV1K5+EKGqYnqYlIbmjL378OTLQbb0GtbGbfSA+vHwthe+C5gmS/yV7H+s8DfAaBlYemXCoyWMwPhfiOdoyRUZ3EuIHrKpxf5xY0dmCiAUz9KkFJjPTG3Tg5la9/pZ9INIeuFLa3DdZIBoS+VNy9r0NaR3UGLYua6iz7BwrLVaCZ/pj75q8RmbRGcmIyR9HQGfTu4ZINL8pvdDx3BikuUuR/wz8BP5MUOG5tAsMAAAAASUVORK5CYII=",[[-33.941641403099,18.3003829808018],[-34.3371886093941,18.5267584324]],1,null,null,"lambda"]},{"method":"addRasterImage","args":["data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAC0AAABfCAYAAAB88E5QAAAFoUlEQVR4nO2a3WscZRTGf2f2I7ubz1ZiSAs1RbGtQrUifoEIpdJGaGsjLUiRgoJ/QL2st8WCF9IbRRRBcuFF0jZabdVGBBEUiuhNWyoFjWBr05C03WyyO7szx4uZ2Z2tTTEz7zS9mAc2mx2yc559ct7znvO8CylSLAlJ4qY6N6oU+7wAhZ3GY1imbzg+8KKGX9sf7NOl/jYqjJL+RLYqQGPiFJSv8Y97hfFtu0yGAAySnj6wQwGGhgpcGq3gnPme87MdfPbDKlMhmjBCWqc/UoBNW0p0loRSj4N7vcrm+1zccw7bDp02miKxSat9SgFKfQ3sukuuoEgGJJ/h7HSGck8Hk4eHjS7GbJw3q31aQaHrfqxSloH+LIObFik8swbZe4CXZGsi1SkWaQDEgmyOji0DrOMqxYOvIv2vC3xqgN4SIaO+0VMZEKH+4cdcnlxkaOKrRJS9FZFzWvLDAorkdojV00GtkjHJK0WKFClSpEiRIkWKFCkMQSvHjbtH/wfxLYSFiRUhHgtaHtOVUjyFzh9TrZxITH3jPoVWjiuqnomDgiqoi3TvNRYrGVO9qXKLNACZHKginXtixTVuqntQEPHUvg29uKmTCGnpHBHchq+wBVbGe6iZNI9vQC4FsQABy/LJegZr3NSAJElDKzWCRYkZpRN1ObU8rliWn9sWUtx1V1zVFClSpLgL0NqXqtWTiffWseqmVk4oIt5tJLT7la/5x3LJIHLvcdsxK7iSycVuiu6ESKSbhC3L7+YItaCC0zfAfK5hkGY7oikddGuqfksR9BXe9Yxdoztvrum/Fcsm3fZvD8hKiJ86xlrQpRCvyxNahCX0+SXZvmjZSkvnHmkjG6SHun5Oa+JKR8ppKe0RXMdP4WC0ksQVDhB93LoTwWCQTQjxNpeg9AX57KeHlOKPVHdCrMFWSi975MI5rYo6k4kmtYFpPORriMUVpr2r9a8TI25gsJW2JBvMrvMtg4QsFUwoHW6YRDzVnXqiZc+MhRAorUA2i2OBJYLrfqfiOnDzKrL6NWOLM7bSUtzdTsapk5EcYlcRp968fGl4pzHpzSSeaqsPcV1wGzT9DnWhu596TVg4NGyEuBHSbqGLVjOtUKu09yKqDD3vIJZQPbI7NnEjeXbTHlNLMhQzPVjVCqCQK3hVpGGD69DsSfyFG8dtMqJ0zirQcG1u2ldZzOGpnMnjCO1bemBKArr4eWTFjVSPgubpyJcQhXJjFrLduDhksP5b+sT/ESNJjJCW3HbRhQlFhG4pAnVqb7/H9SmLwaMj0DMQitjB/i8qzP2UjxzP4Lblq6cOqEt2bReuK8wdGkPP/QyNGqBQr7H90QVytsO+N8Yi6W3Qn9bQk5J96iHWjKxF+tZgvz9KrlqD7k7kwcd4uLfO3OoC6/+8ESmSsV2qzVJoDrrSqtWNGjp1gcaP58k+t5GZDesZvdjBW4/vXzYHc6SDgVdoLbLwNKOuV/UunOXKu78xeHAz1tNHIsU3SPq4tnV8SouwU/eO41zHq93laWTgzZWt0x4CiYOuL7iukPUrhZX1FC/2ovPjkYueMdLelB6SOXSihRsMCXhK5wrt2/wyYfZ0S0P9RyB1m5ET2sY7R+6F9PAPPcON0+2gCvPXcC8eXtltfEmo+inefgBKaRWSLUS+bUIH+se0rXqEv5Hgh4xzcpuYP9Esgd4r2qqLeL9H9UcSTI/wLhO+BucWZ2I5aAmSDlTFey72gl1Bfz/LhkeeJZfZcW9Uj3Z4Smuxl1mrTMOtwfwM1Y1PMFv7O9adEyPtbTYWUp1ndaafRafML47Lyakyl27EC5ug0r7X59TBXqDbtvh1Js+ZqQwNjbf+EyUNIF2vCKp8M1vBUXig1+XEH7lY90ycNHjEL1cs/poXZm2YvR591LrrePKdSX3h6LfptyVTpLgF/wLgRijPgQ+lqQAAAABJRU5ErkJggg==",[[-33.941641403099,18.3003829808018],[-34.3371886093941,18.5267584324]],1,null,null,"A"]},{"method":"addLayersControl","args":[["alpha","gamma","lambda","A"],[],{"collapsed":false,"autoZIndex":true,"position":"topright"}]}],"limits":{"lat":[-34.3371886093941,-33.941641403099],"lng":[18.3003829808018,18.5267584324]}},"evals":[],"jsHooks":[]}</script>